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Abstract 

We calculate the next-to-leading corrections to the jet vertex which is relevant for the 
Mueller-Navelet-jets production in pp collisions and for the forward jet cross section in ep 
collisions. In this first part we present the results of the vertex for an incoming quark. 
Particular emphasis is given to the separation of the collinear divergent part and the 
central region of the produced gluon. 


1 Introduction 


The BFKL Pomeron []^ presents the pertnrbative QCD prediction for the Pomeron, and in 
recent years attempts have been made to verify its relevance for experimental data. Apart 
from the 7 * 7 * total cross section in e’''e“ scattering which is generally considered to be the 
gold-plated BFKL measnrement p|, special jet measurements have been proposed both for 
hadron-hadron colliders (Mueller-Navelet jets |@) and for deep inelastic scattering (forward 
jets [|]). First comparisons of the leading order calculations M with experimental data 
have clearly demonstrated the need of next-to-leading order calculations: both in the 
measurements at LEP and in DIS at HERA the data are below the leading logs (LL) curves, 
while data from the pp collider at TEVATRON are found above the LL estimates. The next- 
to-leading (NLL) corrections to the BFKL kernel ^ lower the theoretical prediction, but 
they are so large that they might even cause serious problems for the stability of the series. 
Various attempts have been made in order to improve the predictivity of the NLL BFKL 

approach. However, so far one has not been able to perform a consistent NLL analysis of the 
data. First of all, a consistent NLL framework for describing not fully inclusive processes, such 
as jet observables, has not been established yet. In addition, even adopting the LL high energy 
factorization formulae at NLL level, a few important pieces of the NLL calculations are still 
missing. 

The three measurements for searching high energy QCD dynamics are illustrated in Fig. ^ 
for a complete NLL analysis one needs, in addition to the NLL calculation of the BFKL kernel, 
the photon impact factor and the jet production vertex. Whereas the former one is currently 
being investigated by two groups |]^,|^, the latter one, so far, has not been calculated. It is 
the purpose of this paper to present a consistent factorization formula for high energy jet cross 
sections at NLL level, and to obtain hrst results for the jet vertex, namely the quark-initiated 
vertex. The case of an incoming gluon will be presented in a forthcoming paper . 
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Figure 1: Schematic diagrams representing the three processes searching for high energy QCD 
dynamics: (a) double DIS or 7 * 7 * scattering; (b) DIS with forward jet; (c) hadron-hadron 
scattering with dijets or Mueller-Navelet jets. 


The theoretical challenge in performing a NLL calculation of such a jet vertex is that it lies 
at the interface of collinear factorization and BFKL dynamics: in Fig. [^c the lower incoming 
parton emerges from the proton and produces a hard jet, thus obeying the rules of DGLAP 
evolution [Q and collinear factorization. Above the vertex, on the other hand, one requires 
a large separation in rapidity between the two outgoing jets (which are assumed to have large 
transverse energies of comparable magnitude). The kinematics between the jets, therefore, 
belongs to large-logs dynamics and is described in terms of the BFKL language. When com¬ 
puting next-to-leading (NLO) corrections^] to the jet production vertex, one expects to find 
collinear divergencies which have to be counted as higher order corrections of the incoming 
parton density; at the same time, part of the NLO corrections will overlap with high energy 
gluon radiation between the two jets which belongs to the leading logs (LL) BFKL approxi¬ 
mation. It is one of the main goals of our analysis to show that both types of contributions can 
successfully be identified and separated from the NLL jet vertex. 

As the central part of our calculation we will compute, to order the high energy limit of 
the cross section of the processes q + q —>■ q + X + jet (Fig. j^), where X may contain one gluon 
or quark. The LL approximation of the order Cl(Q;g log s/sq) has been calculated before [^, 
we will present the NLL (constant in s) term 0{a^). We will show that the cross section can 
be written in a factorized form: there are NLL corrections to the impact factor of the upper 
incoming quark (which have been calculated before 0 ). and for the emission of the gluon 
in the central region we recover the LL BFKL result. Finally, the NLL corrections to the jet 
vertex of the lower incoming quark are what we obtain as new result. Making use of this high 
energy factorizing, we can use our results for the Mueller-Navelet jets (FigjTjc): we can apply 
the NLL results not only to the lower jet vertex but also to the upper one. For the forward jets 
in DIS (Fig. |I|b) we have the NLL corrections for the jet vertex, but we have to wait for the 
completion of the NLL corrections to the photon impact factor. 

An important result of the NLL calculation is the dependence upon the energy scale sq: 
in the leading logs approximation this scale is undetermined and thus introduces a principal 
uncertainty of the theoretical prediction. The NLL calculation determines how the cross section 
changes with a change in sq and removes this uncertainty. 

Moreover, at next-to-leading order, the renormalized parton densities start to play a role, 
reducing drastically the dependence on the collinear factorization scale fip, otherwise maximal 
in all previous LO calculations. 

Our paper will be organized as follows. We begin with a brief outline of our program. In 

^Note the difference between LL (leading log s) and LO (lowest order) and that between NLL (next-to-leading 
logs) and NLO (next-to-leading order = next-to-lowest order = one-loop). 










Secs, ^and ^we present the results of our calculation: first the virtual corrections, then the real 
corrections. Our main emphasis will be on the separation of soft and collinear divergencies at 
the vertex of the lower incoming quark and on the removal of the central region of the produced 
gluon. Combining the real and virtual corrections we obtain in Sec. ^ an analytic expression 
for the jet vertex. We conclude with a brief summary and discussion. 


2 High energy factorization 

2.1 General framework 

The processes that we are going to study are those in which a hadron H strongly interacts 
with parton 6; to be dehnite, we chose 6 to be a quark. In the hnal state a jet J (in the 
forward direction with respect to H) is then identified (see Fig. 0). Our notation uses light 
cone coordinates 

:= ^ ^ , ( 1 ) 

where the light-like vectors pn and pb form the basis of the longitudinal plane: 


Ph = 



s-= (PH + Pbf 


(2a) 


Pb = 



Pi = Ei 


0 Vi 



(2b) 

(2c) 


In the last equation we have introduced a parameterization for the f-th particle in the hnal 
state in terms of the rapidity (in the pn + Pb center of mass frame), of the transverse energy 
Ei = \Pi\ and of the azimuthal unit vector || 



t 


Figure 2: Diagrammatic representation of the high energy process studied in this paper. H is 
the incoming hadron providing parton a with distribution density /; b is the other incoming 
particle which will be taken to be a quark; J denotes the jet produced in the forward direction 
(w.r.t H) and i is the generic label for outgoing particles. 


As usual, our jet consists of a certain number of partons, whose rapidities and azimuthal 
angles are found inside a given (small) region in the {y, 0) plane. The position of the center 
of that region dehnes the rapidity yj and the azimuthal angle (f)j of the jet, its size is related 









to the jet radius R, and the sum of the transverse energies of the particles forming the jet 
constitutes the jet transverse energy Ej. We will keep our jet dehnition rather general; it only 
has to obey a few minimal requirements (see Sec. B- 

We will be interested in the high energy limit of the subprocess a + b^q + i+ jet. In 
addition to the energy s, we need to dehne the momentum transfer 

t-={PH-pjf ■ (3) 

and the jet energy Ej\ 

( qVj Q-yj \ 

71’ 7^’ ^ VsEje-y-^ . (4) 

The high energy (Regge) limit that we are considering is dehned by 

Ej hxed , t hxed , s —*• oo (5) 

which determines a logarithmic growth of the jet rapidity with s: 

l/j~ilogs. (6) 

The use of perturbation theory is justihed because we consider all kinematic scales much greater 
than the QCD scale: 


Ej^ t 3> Aqqq; . (7) 

This condition allows also us to neglect the masses and the Fermi motion of the light partons 
inside H. Equivalently, in terms of the twist expansion, the leading twist contribution can be 
extracted by considering the parton massless and collinear to H. We shall therefore adopt 


Pa = xpH = 



( 8 ) 


It is well known that the Regge limit is dominated by gluon t-channel exchanges |]^ and 
that, in the leading logarithmic approximation, the elastic scattering amplitude and the total 
cross section can be written in a factorizing form: (i) a gluon Green’s function which describes 
the exchanged system, and (ii) impact factors which denote the coupling to the scattering 
partons or particles. NLO corrections that have been computed for the BFKL kernel § 
and for a few impact factors |l6Hl^ support this factorizing form also in next-to-leading order. 
It is one of the goals of this paper to verify this factorization also for the inelastic process 
a + h^a + i+ jet. Thanks to this property we will be able to treat, on the same footing, 
several classes of processes. Two of them are crucial for the study of QCD in the Regge limit: 


dijet, or Mueller-Navelet jets, coming from hadron-hadron collisions where a jet is detected 
in the forward direction of each hadron; in this case particle 6 is a hadron; 


• forward jet, coming from lepton-hadron collisions where a jet is detected in the forward 
direction with respect to the hadron; in this case particle b is identihed with the virtual 
gauge boson, e.g. a photon, emitted by the lepton. 


In the following we assume that a proper dehnition of the jet has been chosen. This choice is 
represented by a function (actually a distribution) Sj which selects the hnal state conhgurations 




contributing to the observable we are interested in. The jet cross section is given by the action 
of Sj on the full exclusive cross section in D = 4 + 2e: spacetime dimensions 

^ oo n 

dcr = ^ y^p7l)^6^{pH +Pb- y^^Pi){\MHb^n\'^) d<h^(pi, • • • ,Pn) (9) 

n=2 i=2 


as follows: 


da da 

dJ dyjdEjd(()j 


daSj. 


( 10 ) 


Here J = {yj, Ej, (pj) collects the jet variables, n is the number of particles in the hnal state, 
M. is the invariant amplitude and d<h is the phase space measure. 

In order to describe perturbatively a hadron-initiated process, we assume — according to 
the parton model — the physical cross section to be given by the corresponding partonic cross 
section dda (computable in perturbation theory) convoluted with the distribution densities fa 
of the partons a inside the hadron H. The partonic distribution functions (PDFs) fa'.aEH 
constitute a non-perturbative input. This approach is justihed provided the infrared singular¬ 
ities stemming from QCD interaction among massless objects can be consistently absorbed in 
a redehnition of the PDFs according to the well known factorization of mass singularities [^. 
Those “renormalized PDFs” will be eventually interpreted as the universal objects measured 
in hadronic collisions and obeying the DGLAP equations. Let us therefore write 

da = V/' dx ddbaix) f^^\x) (11) 


where x = pf/pfj is the longitudinal momentum fraction of the parton a with respect to the 
parent hadron H. In Eq. o we show explicitly that, according to the previous discussion, 
the PDFs must still be considered as “bare” quantities. In conclusion, the jet cross section is 
given by 

^ /dx d&baix) Sj{x)f^^\x) (12) 

a 

and is diagrammatically represented in Fig. ^ 

We proceed by reviewing parton-hadron scattering at lowest order. 


2.2 The Jet vertex at lowest order 

In order to evaluate the jet cross section in the high energy regime 


for parton-hadron 


scattering, thanks to Eq. (JT^), we need only to consider parton-parton scattering. At lowest 
order (LO), the relevant cross section is dominated by one gluon exchange in the t-channel, as 
shown in Fig. Let us dehne the gluon momentum and the partonic center of mass energy 
squared 


k-.= Pi-pa = -zpa+wpb + k± , fc_L = (0,0,fc), (13) 

s:= {Pa+Pbf = XS . 

The partonic cross section is constant in s and is given by 

dd® = hf\k)h^^\k)dk , dk = d^~^‘^'^k , 


(14) 




in terms of the LO partonic impact factors 




h. 


(0) 


(k)-.= N^ 




2i+f 


Qfo 


p?’T{,l-e),/NT^ 


{i = a, h) 


(15) 


where the colour factor Ci'is Cp = — ^)/2Nc for a quark {i = q) and Ca = Nc for a gluon 

{i = g)] the coupling cks in D = 4 + 2e dimensions is dehned in Eq. (^). 

It is evident that the jet can contain only one of the two particles in the hnal state, the 
other moving in the opposite direction. Furthermore, since we are looking for the jet in the 
forward direction with respect to if, the conhguration p 2 = Pj gives a negligible contribution 
to the cross section. This can be easily seen by noticing that the corresponding amplitude 
involves a propagator ~ I/ImI — 1 /s much smaller than that of the pi = pj amplitude 1/t. 
This means that, at lowest order, the jet momentum has to be identihed with pi, so that the 
jet distribution for two-particle hnal states, according to Eq. (0), reads 


sf\puP2]Pa,Pb) = S{yi - yj)6{Ei - Ej)6{(j)i - (j)j) . (16) 

As independent variables in we can adopt = k for the hnal state, and pt/pti = ^ 
the initial state, so that we can dehne 


Sf(k-,x) 


S?\puP2iPa,Pl,) 



(17) 


being kj the transverse momentum of the jet. By substituting Eqs. (H) and (0) in Eq. (|0), 
the LO jet cross section assumes the factorized form 


dcr 

dJ 


E 


dx I dS; )i[°’(fc)/ii“>(fc)5y(fc;i)/™(i) . 


(18) 


Besides the PDF fa and the partonic impact factor hj,, we are left with a term that can be 
interpreted as the LO jet vertex 

q“>(fe,i):=fcl“>(fe)Sf(fe;x). (19) 


The lowest order formula for the jet cross section is therefore 


dcr 

dJ 


[ dx [ dk hP(k)V^°\k,x)fP(x) 


( 20 ) 



Figure 3: Leading diagram at lowest order for parton-parton scattering: the interaction occurs 
via gluon exchange in the t-channel. 











2.3 One-loop analysis: LL approximation and future strategy 

Moving on to higher order, let us hrst address the leading logarithmic approximation, i.e., 
terms of the order Q;”log"'(s/ — t). The resummation of these terms, also referred to as leading 
logarithmic (LL) approximation, was addressed long ago for fully inclusive processes [^, and it 
has also been applied to dijets p[ and forward jets 0. It is instructive to review briefly how the 
logarithmic enhanced terms arise at one-loop, and to study the structure of the singularities 
and their connection with the various kinematic regions. 

b .^2 

ki 
ki 



Figure 4: Labelling the two-to-three parton scattering process. Note the dehnition of the 
transferred momenta k and k'. 

We consider hrst the real corrections to, say, quark-quark scattering, which involves the 
emission of an additional gluon of momentum as shown in Fig. The structure of the hnal 
states giving the leading contributions corresponds to the so called multi-Regge kinematics 
(MRK) 

2/1 > 2/3 > 2/2 , .Fi ~ .Fs ~ -^2 , ( 21 ) 

where the rapidity 1/3 of the emitted gluon is strongly ordered between the rapidities of the 
scattered partons yi and 1/2 (i.e., the gluon is emitted in the central region), while the magnitudes 
of the transverse momenta are of the same order. 



Figure 5: Feynman diagrams contributing to the qq qqg process at leading logs order. 


All diagrams shown in Fig. contribute at LL, and the resulting differential partonic cross 
section reads 

1 




-dkdk' — 

TT 7r,(fc - k'f Z ’ 


( 22 ) 


where vr^ is given in Eq. ([^). We have introduced the momenta transferred by the a and b 
partons 

k:=pb-p 2 =-wpa + wpb + k± , k± = {0,0,k) (23a) 

k'-.= Pi - pa =-zpa +zpb +k'^ , = (0, 0, fc') . (23b) 

The transverse energies introduced in Eq. (^ correspond to Ei = \k'\,E 2 = |fc|, E 3 = \k — k'\. 
In the MRK (pi]), so that the outgoing parton 1 carries most of the 

plus-component of the parent parton a. By using Eq. (^ and comparing with Eq. (H) we get 

Xy/S 


2/1 - log ■ 


El 


(24) 
















Since Ei = |fc^|, k' characterizes completely the momentum pi. The range of the momentum 
fraction z is approximately given by pi/s < z < l — k''^/s. The boundary values of ^ correspond 
to the gluon in the fragmentation region of parton b and a respectively, where Eq. no longer 
holds. Nevertheless, it shows that, upon ^-integration, a logarithmic factor log(s/so) appears: 


dd 


(l,real) 

ba 



(l,real) 

ba 






1 

TT 7rs{k — k'Y 


;^{0,real)(yj,,y^/) log A + cOUSt 
■^0 


dfc dk' 


(25) 

(26) 


The coefficient of the logs term jg ^i^g leading log BFKL kernel. The 

scale of the energy sq is a parameter of the order of the transverse momenta squared of the gluons 
(so pI = |fc — ~ ~ k''^). Its value is not determined in the LL approximation — a 

change of sq would affect only the constant piece in Eq. (^5[) — , but will play a central role when 
going to the next-to-leading level of accuracy. In particular, that part of the differential cross 
section in which the emitted gluon lies outside the central region contributes to the constant 
term (i.e., without a logs enhancement). 

Let US now consider the virtual corrections. In covariant gauges, the diagrams involving 
two gluon exchanges (see Figs. ^,b) give contributions to the amplitude which increase loga¬ 
rithmically with the energy, whereas other diagrams (shown in Figs. ^,d) have no logarithmic 
enhancement. Note, however, that the latter diagrams contain the ultra-violet (UV) singulari¬ 
ties that provide the renormalization of the coupling. 






Figure 6: Feynman diagrams contributing to the one-loop virtual corrections to qq —> qq elastic 
scattering in LL approximation (a,b) and in NLL approximation (c,d,...). 


The result of the virtual correction to the partonic cross section can be presented as follows: 




2u^^\k) log + const 

k 


dk 


(27) 


The coefficient in front of the log s is the one-loop Regge-gluon trajectory, because — 
according to Regge theory — the amplitude for ab —>• ab elastic scattering mediated by color 
octet exchange is described by the exchange of a reggeized gluon and takes the form: 


Mba ~ 



= 1 +aai.j6)(t) log^ H-. 


(28) 


By taking the square modulus and adjusting the normalization through the LO expression (|T^, 
it is straightforward to match Eq. (0)- This contribution to the cross section dehnes the virtual 
part of the LL BFKL kernel: 


A'l“''’'">(fe, k') - 2ui<-'\k}S(k - k') . 


(29) 
















Together with Eqs. and (^) we obtain the fnll one-loop partonic cross section 




(fc, k') log-h const 

So 


dfcdfc'. 


Here collects the real and virtnal coefficients of the logarithmic terms: 

k') := k') + k') . 


(30) 


(31) 


and dehnes the LL BFKL kernel; the term const contains contribntions that have no large 
energy logarithm. 

Let ns nse these resnlts to dehne the jet cross section. The jet we want to observe inclndes 
the ontgoing particles carrying the largest valnes of rapidity lying within a specihed (small) 

range R. However, the strong rapidity order in Eq. (plf ) allows only pi to enter the jet, so that 

( 3 ) 

the jet distribntion for three particles in the hnal state S) rednces to the lowest order one 
sf (Eq. (|I^) ) which, combined with the impact factor in Eq. (^OD, reprodnces the lowest 
order jet vertex (Eq. O). Convolnting with the PDF we obtain the following factorized 
expression for the one-loop LL jet cross section (s = xs)\ 

^ j dx j dkdk' hf\k)K^^\k,k')\og^V^^\k',x)f^\x) . (32) 


It is the pnrpose of this paper to compnte, in this order of cts, the constant (i.e., non- 
logarithmic) corrections to onr jet cross section formnla. We generalize the resnlts (^) and ([3^ ) 
and make the following ansatz: 


^ = q ; s ^ j dx j dkdk' hb{k)G{xs,k,k')Va{k',x)fa{x) , 


G{xs, k, k') := 6{k — k') + (k, k') log — , 

So 


(33a) 

(33b) 


According to the previons remarks, we suppose that the inclnsion of the one-loop constant 
terms just provides perturbative corrections to the quark impact factor, to the jet vertex, and 
to the PDF as follows: 


h = 

V = 

f = . 


(34a) 

(34b) 

(34c) 


Equivalently, our ansatz corresponds to the following structure for the one-loop cross section: 
dad) 


dJ 


— Q^s 


j dx J hf\k)K<'»{k,k')log^Vi°\k',x)ff\x) 


(35) 


;4‘>(fe)l4<“>(fe,l)/i“>(l) + ;4"'(fe)q“>(fc, + hf(k)V!l'\k, x)fi»\x) \ , 


which is obtained simply by expanding Eq. ( |33a)) up to relative order ag. 

In Eq. (^) the Born approximations (marked by the superscript (0)) have been listed 
in (0), d^), and (|g) respectively. For the hrst order correction to the partonic impact factor, 
hd), which appears on the second line, we can use the known expression of Ref. |^, ^, and 












for the correction to the PDF, we have the usual LO Altarelli-Parisi splitting functions 
(in the MS scheme): 


27re: 

1 




27r£ \ fi 


/i" 


E ® /P ■ 


(36) 


Finally, the correction term is what we want to compute in this paper. 

Eqs. (^) and (|^) constitute a highly non trivial ansatz, which will be shown to depend 
upon a careful separation of singular and hnite pieces. Our main task will consist to identify 
the collinear singularities (|36D which go into the renormalization of the parton densities, to 
check that the other infrared singularities cancel out when adding virtual and real corrections. 


and, hnally, to separate the terms proportional to logs which belong into the first line of (|35[) . 
At the end, a hnite term remains, which eventually can be interpreted as one-loop correction 
to the jet vertex, Va^\ 

In the rest of this paper we will compute this jet vertex correction, Va^\ In this paper we 
will concentrate on the case of incoming quark {a = q). The case of an incoming gluon (a = g) 
will be dealt with in a subsequent paper |l^ . 


3 Virtual corrections 

In the following we will develop the one-loop analysis of the quark-initiated jet production 
process. We adopt dimensional regularization in D = 4 -|- 2e dimensions and dehne, according 
to the MS scheme, the dimensionless coupling as as a function of the dimensionful bare coupling 
g and of the renormalization scale jj, as follows: 


tts = : = 


9V'r(i-£) 


(37) 


(47r)^+^ 

We begin by collecting the virtual corrections. Some of the diagrams are shown in Fig. ^ 
Discarding all terms suppressed by powers of s, the one-loop parton-parton cross section can 
be derived from Ref. and reads 


ddfea = tts h^^\k)h^°\k) 


2a;d)(fc)log^ + n 6 (fc) 
k 


na(fc) 


dk 


(38) 


The hrst term has already been introduced in Sec. |2^: it represents the LL contribution to the 
virtual corrections. In particular the coefficient of logs, namely 2tnd), constitutes the virtual 
part of the leading kernel of Eq. (pTI) and is just twice the one-loop Regge-gluon trajectory 


u 


( 1 ) 


{k) = 


Ca 1 F2(1 + £) /k 


(39) 


71 2e F(1 -I- 2e) 

It shows an £-pole due to a soft singularity which compensates the corresponding one of the 
real part of the kernel. 

The non logarithmic terms in Eq. (^) represent the NLL contribution to the virtual correc¬ 
tions and are expressed in terms of the virtual corrections to the impact factor 11. The virtual 
corrections to the quark impact factor read: 


n,(fc) = 


llN^ - 2Nf 1 
£ 


+ 


85 TT^ 

^ ^ T 


TT 


18 TT 


-h4- 

26 


TT 

”6 


Cf 

TT 


/i 2 


(40) 














In the above expression we have singled out terms of different physical origin. The hrst term is 
proportional to the /3-function coefficient bo = — 2Nf)/12n. It multiplies the ultraviolet 

(UV) pole providing the MS renormalization of the coupling 


as(fc^) := 




£ \/i 


2 \ £■ 


(41) 


In fact, at LO the partonic cross section (|T^ is simply the product of two bare partonic impact 
factors 


J - (0) 


(42) 


where we have explicitly shown only the dependence on Adding the UV divergent term 
of Eq. (|38D stemming from 11^ renormalizes the coupling inside hq\ 


dd 

dk 


( 0 ) 

bq 


+ 


dd 

dfc 


(1,UV) 

bq 


Uo 




1 - a!“> 


= . (43) 


The same UV pole can be found in Ilf,, and it provides the running of the coupling of 

The virtual contribution to the jet cross section is easily obtained by substituting in Eq. ([l^) 
the expression (0) for the jet distribution and Eq. (|38|) for the partonic cross section. By 
combining the jet distribution with the quark q impact factor at LO we reproduce the LO jet 
vertex (^) and, after renormalization, we end up with 


da 


dJ 

where 


- Qo 


dx I dk hf\k) 


2j^\k)\og^ + n,ik)+Uqik) 

k 


(44) 


n,(fc):= n,(fe) - (- 6 „/£) 


85 


71 


c, 


5 Nf 


71 


^ LT I - Ts -— t )-^0 log — 

36 d J 71 18 71 \e^ 2e 6 / tt 


C, 


k 


2 n 


e 

/i 2 


(45) 


Any occurrence of ag in Eq. (|4^) and in all other coming formulae is to be understood as as(h^)- 
The “reduced” quark impact factor virtual correction (]45|) shows double and single poles in 
£. These poles are of IR origin and are due to both soft and collinear singularities. Partly they 
will cancel against the corresponding singularities of the real emission corrections, leaving a 
simple pole that will be absorbed in the redefinition of the PDFs. This will be shown in Sec. ^ 


4 Real corrections 

When calculating the real emission corrections to the one-loop jet cross section, our main 
concern will be the correct treatment of the IR singularities (keeping in mind that we are 
considering a partially exclusive process). Infrared singularities are contained in the upper 
quark impact factor, the real part of the BFKL kernel, and in the lower jet vertex. Some of the 
latter ones contribute to the renormalization of the incoming parton density. A priori it is not 
evident that the overlap of the various regions of the phase space responsible for the divergencies 
can be disentangled in such a way that they reproduce all the expected singularities (and not 
more). We will show that this is actually the case. 















4.1 Jet definition 


We begin with a brief review of the jet dehnition. We follow the arguments given in 
and we wish to keep our distribution functions S) as general as possible. In massless QCD 
two kinds of IR singularities exist: (i) soft singularities which arise when a gluon is emitted 
with vanishing momentum; (ii) collinear singularities which arise when two interacting partons 
are emitted collinearly. In order to dehne infrared hnite jet cross sections we have to require 
hnite limits whenever momenta in the hnal state belong to either (i) or (ii). Given a set of 
functions ,Pn]Pa,Pb) (where (pa-iPb) denote the momenta of the initial state), we 

have to require that a state (pi, • ■ ■ ,Pj, - ■ ■ ,Pn) with a soft particle pj ^ 0 be indistinguishable 
from the n — 1-particle state (pi, • • • ,Pn): 

\imSf\--- ,Pj,--- ]Pa,Pb) = ■,Pa,Pb) (46a) 

Pj^O 

with Pj being dropped in the RHS. In the same way, the n-particle state (pi, • • ■ , Pi, Pi+i, • ■ ■ , Pn) 
with two collinear particles, e.g. pi || Pi+i, cannot be distinguished from the n — 1-particle hnal 
state (pi, • ■ ■ ,Pi + Pi+i, ■ • ■ ,Prt)- The jet function must then fulhll 

,ap,bp,--- ]Pa,Pb) = ,(a + 6)p, ;Pa,Pb) , (a, 6 > 0) . (46b) 

When an outgoing particle is collinear to an incoming one, say a, a similar relation (which can 
be inferred from Eq. ( [16b|) by invoking crossing symmetry) holds: 

S'j”^(-■ ■ ,apa, ■ ■ ■ ;Pa,Pb) = ; (1 - a)pa,Pb) , (0 < a < 1) . (46c) 


The last equation expresses the property of factorizability of initial state collinear singularities. 
Eqs. ([46b|) and ( |46c| ) should be understood as smooth limits for momenta approaching the 
collinear conhguration. 

In our case real emission involves three partons in the hnal states, one gluon in addition 
to the incoming quarks a and b: abg. As indicated in Fig. we label the outgoing partons 
in our process ab 123 by 1 = (quark a), 2 = (quark b), and 3 = (gluon g). Let us start 
by listing the possible IR singular conhgurations. Only the emission of gluon 3 with vanishing 
momentum gives rise to soft singularities. Collinear singularities arise in collinear emissions 
of partons that couple directly to each other, i.e., have a common vertex. This is the case 
for pairs of gluons, for quarks and gluons, and for identical incoming and outgoing quarks (or 
quarks and antiquarks). Therefore, the list of all possible collinear singular conhgurations, in 
an obvious notation, reads as follows: 


a II 1 , a II 3 , 1 II 3 , (47a) 

6 II 2 , 6 II 3 , 2 II 3 . (47b) 

It is important to note that, in the kinematic regime we are considering, conhgurations in which 
quark 2 is emitted outside the fragmentation region of quark b are strongly suppressed, i.e., 
quark 2 never belongs to the jet produced in the forward direction of quark a. To see this 
we note that as long as quark b is dehected at small angles by gluon exchange, the amplitude 
contains a 1/f factor coming from the gluon propagator. Suppose, on the other hand, quark 
b dehected with an angle large enough to enter the central region or even the fragmentation 
region of quark a (which includes the jet region). If the scattering happens by gluon exchange, 
the propagator of the latter is of order 1 /|m| — 1 /s, providing a suppression factor ~ t/s 
compared to the small angle case. There is also the possibility of quark exchange that involves 







a 1 /t propagator, but in this case the spin of the quark introduces a new suppression factor 
~ —t/s compared to the gluon exchange. Therefore, we can safely neglect the conhgurations 

in which quark 2 enters the jet, and only particles 1 and 3 play a role in building up the jet. 
To become more specihc, we change the argument structure of the jet distribution functions 
and introduce, as independent variables, Pi,P 3 ,P 3 


= Sf\k',k- k\xz-x) := ^j^^(pi,P 2 ,P 3 ;Pa,Pb) • 

Ph Ph 

:( 3 ) 


(48) 


The soft IR constraint (^ 6 |a) of Sj applies only to soft gluon emission and, since quark 2 does 
not participate in the jet, the collinear conditions (^b,c) apply only to the conhgurations listed 
in (|47a)) . The corresponding relations for read: 


3 soft : 

‘5®(P,0,0;a;) 


(49a) 

111 3: 

Sf\ap,bp,^;x) 

= + b)p-, x) 

(49b) 

a 1 1 : 


= sf\p;0 

(49c) 

a 1 3 : 

o' 

= sf\p-,x-0 ■ 

(49d) 


In the following sections these relations will be used when extracting the divergencies of the 
real emission. It is crucial that the singular contributions generated by the real corrections are 
proportional to the LO cross section: only in this case cancellations with the virtual singularities 
can occur, and the factorization of the collinear singularities into the PDFs can be performed 
consistently. Therefore, the reduction in the IR singular conhgurations contained 

in Eqs. (^91) is a necessary prerequisite. 


4.2 Phase space splitting and the master formula 

Before embarking in the analysis of the one-loop real corrections, let us divide the phase space 
of the outgoing gluon and present the ’master formula’ that we are going to make use of. We 
have already pointed out that the log s term arises from the conhgurations with gluon 3 emitted 
in the central region, while the gluon in the fragmentation region of b should mainly contribute 
to the impact factor correction whereas the gluon in the fragmentation region of q should 
provide the jet vertex and the PDF corrections fq^\ It is easy to dehne a rapidity cut that 
separates the two fragmentation regions: we perform a boost into the positive z-direction which 
takes us into the partonic center of mass system (PCMF). It shifts rapidities while transverse 
energies and azimuthal angles are preserved. Consequently, a four momentum is transformed 
into p'^ with 

P'^ = = (e^^p+,e“^V,P) , A?/ = ^log- , (50) 

2 X 

and rapidity is shifted according to y' = y + /S.y. In the PCMF, we dehne the cut as the rapidity 
center = 0 which corresponds to ycnt = the Sudakov parameterization (p3D, 

the form of the cut is very simple: 


2/3 = 2/cut = 0 


^cut - -^cut 



(51) 


Correspondingly, our rapidity phase space is divided into the “upper half region” (negative 
rapidity: 2/3 < 0 or tc > which contains the fragmentation region of quark b and half of 






the central region, and the “lower half region” (positive rapidity: 1/3 > 0 or z > which 
contains the other half of the central region and the fragmentation of parton q and contributes 
to the jet vertex. 

Finally, we need the partonic differential cross section for the bq —>■ 123 process. They have 
been computed in Refs. In the high energy regime, where we neglect terms suppressed 

by powers of s, the form of the partonic differential cross section turns out to be quite simple 
when restricted to one of the two halves of the phase space, j/3 < 0 or j/3 > 0. For the “lower 
half region” 2/3 > 0 the cross section can be cast into the general form 

daba^fin = hf\k)Ffinik, k\ z)h^^\k') dk dk' dz , {z > z^ut) , (52) 


where the function F depends on the particular hnal state. In particular, for quark-quark 
scattering, we have 

F,,g(fc, k',z) = 1 [CrZ^k'^F CA{l-z)q-{q- zk)] , (53) 

Vgq{z, e) = , TTg = 7r^+^F(l - , (54) 

where q = k — k' is the gluon transverse momentum and Vgq{z,e) is — apart from a missing 
Cp factor — the real part of the q ^ g splitting function in 4 -|- 2e dimensions. In the “upper 
half region” 7/3 < 0, the same relation ( |52[) holds with the replacements 

k —k' , k' —>• —k , q q , z ^ w , (55) 


except for the impact factor which retain their form. These ’master formulae’ (|^ - (|55D will 
be used in the following in order to hnd the real corrections to the quark-initiated jet vertex. 


4.3 Real corrections to the upper quark impact factor 


We begin by computing the contribution to the jet cross section given by the “upper half region” 
j/3 < 0 which gives rise to the NLO impact factor of the upper quark. The starting formula 
is derived from Eq. ([T^) , using Eq. (H) for the jet distribution (with 2; = E^/wxs to good 
accuracy) and Eqs. (|52D and (^) (with the replacements (^)) for the partonic cross section: 

c1(t(2'3<0) as /■ ,. ,, , , mu, X wn^ ,, ^^Vgq{w,e) 1 


dJ 




dkdk' 


(56) 


' R’cut 


TTg q'^{q + wky 

rl Tp2 

+ CA^i-— w)q ■ iq + wk')\ / dxS^?{k\q, — —x]x)fj;^\x). 

Jo 

Since 2/3 < 2/3 < 0, the gluon is emitted very far from the jet region ?/ ~ j/j ± R, for which 
2/j ~ log a/s > R 1 (cf. Eq. d^)). Therefore, only quark 2 can enter the jet, so that in this 
half of the phase space 


Sf{k',q, 


x;x) = Sf^{k',q, —;x)= Sf\k';x) , (w > Wcut) 


(57) 


wxs ws 

It is now possible to factor out from the w and k integrals the q impact factor and the jet 
distribution, which, according to Eq. ( 0 ). reproduce the LO jet vertex. We obtain: 
j (i/^<0) r f 

= as dx dk' U{k',x)V„^°\k',x)jf\x), (58) 


dJ 
U{k',x): = 


J\fCf 

271 


dk 

TT^kJ 


dta 


' R’cut 


w 


[1 + (1 - wf + ew^] 


Cpw'^k'^ + Ca (1 — w)q ■ {q -|- wk') 
q^{q + wk'Y 


(59) 











The computation of the U integral can be done following the calculation of Ref. [^. We 
repeat here the main steps. Let us consider separately the two terms involving different colour 
constants. 

The integrand of the Cp term of U is regular — actually vanishes — for tc —>■ 0 at fixed k, 
so that the lower bound Wcut can be set equal to zero, introducing a negligible error of order 
/ s. Changing the transverse integration variable k q = k — k' yields 


Uc, 


27r Jo 


wY + ewY 


dq 1 

Tie Q^{q + wk'Y 


(60) 


Until now, we have not mentioned any upper limit of the transverse momentum integrations. 
We know that the transverse momentum squared of the outgoing particles are kinematically 
limited to some value of the order of s. For our purposes, however, these upper limits are 
not important, because the transverse momentum integrals converge in the ultraviolet region 
fc —>■ oo, and extending the upper limit of integration to infinity causes negligible errors of order 
1 /s: 


r ^ 

iqi~\/s q 


oc 




1 

s 


The transverse integral in Eq. (|60|) can then be easily performed. 


f ^ 1 

J Tie q‘^{q + wk'Y 


w2e-2 /fc' 2 y r 2 (£) 
fc'2 [jE) r(2£) 


w2e-2 /fc/ 2 Y n 

fc '2 \ /i 2 / e 


6 


e + 0{e‘^) 


(61) 


(62) 


The £ pole reflects the two collinear singularities at q = 0 and q + wk' = 0 which correspond 
to 6 II 3 and 2 || 3, respectively. At this point also the tc-integration can be easily performed. 
Because of the pre-factor in Eq. (|B2[), the tc-integration turns out to be divergent^ at 

w = 0. This is just the soft singularity expected for > 0. The final result for Ucp is 


Ucp 


—hf(k’) 

TT 



1 3 ^ TT" 

~ ~ + 4 —- 

2e 6 


(63) 


up to terms 0(e). For the impatient reader we note that the whole Ucp cancels out when added 
to the virtual corrections of Sec. 

^The reader could object that the w = 0 singularity is an artifact of having set Wcut = 0 which corresponds 
to having neglected an infinite contribution. However this procedure is consistent in that the logic of this 
calculation is first to perform the Regge limit (^ at e 0 and only as a last step, after cancellation of the 
divergencies, to perform the physical e ^ 0 limit. 



Figure 7: Structure of the Feynman diagrams contributing to the qb —> qbg process in NLL 
approximation when the outgoing gluon is emitted with negative rapidity t /3 < 0. The blob 
represents all possible bg* —> bg QCD sub-diagrams. 













The Ca term of U requires special care, because it contributes to the high energy leading 
logs piece of the cross section. In Sec. (Eq. we have presented the leading partonic 

differential cross section, which is nothing but the partonic cross section (^) with Eq. (|53|) 
evaluated in the z —>■ 0 limit (the central region). The same consideration holds also here after 
the replacements (^). In particular we identify the leading term of f/^A with 

UCaCf " ^ 


uh\ = 


TT 


dfc 1 
^ {k-k'y 


dta 


k, k') 


= 2ti I dk 


dtc 


w 

'd(tc, k, k!) 
w 


(64) 


where is dehned in Eq. (^). The function ■d in the above equation signals that the 

functional form of the integrand, accurate in the central region, has to break down somewhere 
in the fragmentation region of quark b. According to the analysis of Ref. ||^, the emission 
probability in the splitting q q'g is dynamically suppressed when the emission angle of the 
gluon g is smaller than that of the quark q'. In the present case, by comparing the ratios of 
transverse to longitudinal components of particles 2 and 3, one expects the active phase space 
to be 


> 




w < 




W 


W 


E/2 + Eq 


( 66 ) 


This has led the author of Ref. to propose a leading term of the form of Eq. (pdD with 
'd(w, fc, k') := 0(^3 — WE 2 ), which matches Eq. (p^) in the low-tc region. With this choice. 




dk h^^\k)K^°’'^‘^^^\k,k')\og 


r(fc, k') 


( 66 ) 


r(fc, k') = max(E 2 , E 3 


The remaining part of Uqa^ i-e., Uca~Uq^, is constant in s, and hence it is a NLL contributions; 
it is given by an integral that, at hxed transverse momenta, is now hnite for tCcut —> 0. The 
result is 





3 TT^ 1 

1 

1 

_“4£ ~Y~4_ 


(67) 


The e-pole stems from the transverse fc-integration in the neighbourhood of the singularity at 
k = 0. 

The complete contribution of the “upper half region” to the jet cross section can be conve¬ 
niently presented if combined with the virtual correction contribution of Eq. (^4|) coming from 
the b impact factor correction 11 ;,: 

ci(T(*'3<0) dcT 


dJ 


— Qfc 


Q!o 


dJ 

dxj 

dx I 


nt 


dA: dk' k') log 

dkhi'\k)Vf>(k,x)ff>(x), 


xs 


( 68 b) 


In addition to a LL part, this formula reproduces the hrst constant term of Eq. (|35|) (of course, 
only the term a = q), namely the full one-loop impact factor correction of the upper quark b: 




C(fc) = 


TT 


3 e 

“4 + 4 


5W 


1 67 TT^ 

e ^ ^ ~ 12 ~ 18Ca 


7 


/i 2 


- bo log ■ 


k^ 


(69) 


























The quark impact factor has been calculated in also in Ref. |]T^, but there a different dehni 


tion has been used. In order to explain the relation between the two approaches, some general 
remarks on the dehnition of impact factors and energy scales might be in place. It is known that 
processes involving coloured incoming particles are affected by collinear singularities that lead 
to divergent cross sections. These singularities depend on the type of the incoming particles, 
and it is therefore natural to associate them with the process dependent impact factors. In our 
paper, we follow Ref. MM and we require the partonic impact factors to include singularities 
of collinear origin only. Such a prescription may sound somewhat academic, because partonic 
impact factors have no phenomenological application and interpretation. However, in view of 
the jet vertex which in the present paper represents our main goal, this requirement is very 
natural. 

Namely, at the lower end of the diagram, where the coupling of the reggeized gluon to the 
incoming parton is described by the jet vertex (more precisely, by the convolution Va® fa which 
we may call “jet impact factor”) it is mandatory that its singularity structure matches the one 
required by collinear factorization. Only in this case /„ can be identified as the usual parton 
density with one-loop corrections (|36D , and only with this prescription the remaining jet vertex 
Va is hnite. As we shall see in the next section, the jet distribution function helps to disentangle 
this structure from the leading log part in a very natural way, but the basic requirement is the 
matching of the singularity structure with the collinear singularities. 

In our approach we therefore insist on the collinear properties of any impact factor. It is 
apparent from Eq. (|69|) that has a simple pole which is the one expected from the & || 2 
collinear divergence — more precisely, the residue of the pole reproduces the non singular part 
of the 7 gq anomalous dimension. This is a consequence of our choice for the LL subtraction 
needed to dehne the impact factor: the angular ordering prescription takes proper account of 
the whole neighbourhood of the & || 2 collinear region, and it avoids spurious singularities which 
potentially could be present due to the {k — k')"^ denominator in Eq. (| 6 ^. Nevertheless, even 
with this requirement on the singularity structure there is still some freedom in the dehnition 
of the impact factor — essentially a factorization-scheme arbitrariness — which corresponds 
to changes in its hnite (in e) part. It is possible to show that the “good” collinear properties 
of the impact factors are preserved, provided the angular ordering prescription in the leading 
term is fulhlled in the soft (small w, E^) region, while the details of the subtraction at hnite w 
or E 3 ahect only the hnite part. 

Returning to the evaluation of leading term, one might think of two diherent strategies. 
First, we could consider the whole term as part of the kernel; this amounts to perform 
the w-integral in Eq. (|5^ and to obtain the LL kernel, multiplied by a logarithm of the energy. 
This corresponds to the hrst line of Eq. ([ 6 ^), where the energy scale r (the denominator in 
the argument of the log), being a function of k and k', is determined by the particular LL 
subtraction, i.e., by d. If we adopt, for instance, the angular ordering prescription in the whole 
range of w, i.e., 'd(tc, fc, k') := 6(^3 — w{E 2 + E^)), the energy scale is r = E 2 -|- E 3 . In this 
case the new impact factor dihers from Eq. (^) only by a hnite piece. 

Alternatively , we might chose a particular energy scale, say f, and then decompose into 
a leading term containing the \og{^/afs/f), plus a next-to-leading term containing the log(r/f). 
The latter term can be expressed, at relative order O^af), as a multiplicative operator factor: 


-I- ash^^^)(l + log 


HP{k,k'):= K^°'>{k,k')\og 


r{k, k') 


(70) 

(71) 


r(fc, k') 

As long as we have chosen the energy scale f in such a way that in the limit k k' [k — k' k) 







it reduces to f —|fc'| ~ |fc — k'\ = {f ^ |fc| = E2), the factor (1 + cusHl) can be safely 
embodied in the impact factor term without spoiling its collinear properties, i.e., it changes 
only by a hnite (in e) amount. 

On the other side, for Regge-motivated scales of the energy like the one proposed in Ref. 0], 
the inclusion of the term in the impact factors^ leads to a different infrared behaviour, giving 
rise even to double poles. In this case, in our opinion, it will be useful to separate the factor 
(1 + asHi) from the impact factor and eventually to include it into the dehnition of the NLO 
BFKL kernel. 


4.4 Real corrections to the jet vertex 

In this section we consider the real corrections in the “lower half region” 7/3 > 0, i.e., the 
corrections to the jet vertex. The starting formula is derived from Eq. (^), using Eq. (^) for 
the jet distribution and Eqs. (|5^) and (^) for the partonic cross section: 


d(jfe3>o) a 
dJ 




X 


[Cpz^k''^+ Ca( 1 - z)q ■ (q - zk)] f dx Sf \k\ q,xz\x) , 


(72) 


where we have used Eq. (^) in order to write explicitly the kinematic dependence of the q 
impact factor. 



Figure 8 : Structure of the Feynman diagrams contributing to the qb —> qbg process in NLL 
approximation when the outgoing gluon is emitted with positive rapidity y'^ > 0. The blob 
represents all possible qg* qg QCD sub-diagrams. 

In this phase space region, the quark-quark {b, 2} subsystem is kinematically well separated 
from the quark-quark-gluon {q, 1, 3} system, and they are connected only via the exchange of 
a gluon of momentum k (Fig. |^). Therefore, at hxed k, the dynamics of the two subsystems 
are independent of each other. Actually, because of the fc-factorization only the transverse 
momentum k and the longitudinal component w need to be hxed in order to separate the 
two systems, because the {b, 2 } system does not contain hnal state emissions at large sub¬ 
energies, and its dependence on the longitudinal component iv is very weak and can be neglected. 
Moreover, the value of w is constrained by the mass-shell condition of quark 2, so that only k is 
the relevant variable between the two subsystems. We therefore expect a “pure” fc-factorization, 
where the { 6 , 2 } coupling with the exchanged gluon is described simply by the h^\k) impact 
factor. In order to hnd the expression for the {q, 1, 3} system and its coupling to the gluon 
(and to the hadron H), we hx the transverse momentum k, remove the b impact factor and 
study the remaining part of the jet cross section Eq. (ff^). 

^In Ref. lH^, one can find the expression for the impact factor at the Regge-scale f = |fe'| = £'2 and the 
equivalence with our expression -t asH^[^) is proven. 





First of all, it is important to understand the structure of the singularities in Eq. ( 0 - 
The integrand contains three singular points in the fchintegration, namely the zeroes of the 
denominator k' = 0, q = 0, and q — zk = 0. These points correspond to the collinear 
conhgurations g || 1, g || 3 and 1 || 3, respectively. Moreover, there is a potential soft singularity 
hidden in the 1/z pole of the “splitting function” Vgq. The numerators in the square brackets 
of Eq. (^) soften some of those singularities, but this happens differently for the Cp and the 
Ca parts. Therefore, we consider these two terms separately. 


4.4.1 Cp term 


The Cp term, owing to the k!’^ factor in the numerator, has no g || 1 collinear singularity. Due 
to the factor z^ in the numerator, the z-integrand is no longer singular at z = 0 and, along the 
same line of arguments as given in Sec. |4.3|, we can shift the value of Zcut 0. Dropping the 
( 2/3 > 0) label, the Cp part of Eq. (El) is given by 


da 

dJ 


Cf 


= a, 


^ dk hf\k) [ dx X 
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'0 


»1 


d^; z‘^Vgq{z,e) 


dk' 


1 


Sf\k', q, xz] x) 


q^iq - zkf 


(73) 


It is convenient to rescale the gluon transverse momentum by setting q = . zl, and to use I as 
integration variable by substituting k' = k — zl, so that q — zk = z{l — k). Next we perform a 
simple fraction decomposition in order to separate the initial {i) state (g || 3 Z = 0) and 

hnal (/) state (1 || 3 I — k = 0) collinear singularities: 


i\i-kY c + {i-ky 




(Z - ky 


(74) 


Beginning with the final state (/) collinear singularity, in terms of the new variables the 
Cp contribution to the jet cross section can be rewritten in the form 


da^ 

dJ 


:= 


dz 

yl-2s 


dZ da^’"""^* da^’™^^ da^’®”**" 

— ky ’ dJ dJ dJ 


(75) 


which is particularly suitable for the analytic extraction of the divergencies: the RHS contains 
three pieces, 1) the soft divergence, 2) the collinear divergence, and 3) a hnite part. The explicit 
expression of the integrand I in Eq. (E^) is 


I{z,l)-.= ^zVgq{z,e) j dk hf\k) 


MC. 


dxSf\k — zl,zl,xz-,x)fj;^\x) . (76) 


z^ + (z - ky Jo 


The soft term in Eq. ([75|) is dehned by evaluating the integrand in the soft limit 2 ; —> 0. 
In this limit, the jet distribution can be simplihed by means of Eq. (^9^, and leads to the 


constraint k^ = Ej. One obtains 
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dk / dxhf\k)Vj^\k,x)f^^\x), 





















where, in the final result, we have collected some factors in such a way that they reproduce 
the LO jet vertex (|^. It can easily be seen that we have recovered the LO structure of 
the factorization formula. The additional divergent factor exhibits single as well as double 
poles, because our definition of the soft part includes also the region where collinear and soft 
singularities merge. 

The pure collinear singularity can be isolated by evaluating the integrand ( 0 ) in the collinear 
limit I = k, after having subtracted the soft term [I = k, z = 0). The resulting expression is 
clearly regular in the soft limit {z 0) and therefore contains a simple collinear pole. An UV 
cutoff A is introduced since the residue at the collinear limit is no more integrable in the UV 
region. Thanks to Eq. (' |49b|) the jet distribution simplifies to Sf'^: 




dJ 


:=as 


dz 


dl 


J — ky 

UC 
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The remaining part is regular in the £ 


dx [ dk h® (fc)l/^(fc, x) A^(a;) . 
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0 limit and defines the finite term: 


d^hfinite 


dJ 


:= as 


'0 


^ d^ 


dl 


7r(Z — k)'- 


[Hz, l) - /(o, l) - (Hz, k) - HO, k))B(A^ -d- kf-)] 


(79) 


This result (with explicit expressions for I{z, k) etc.) will be used in our final formula. 

Next we consider the term (i) with the initial state collinear singularity. We can write, in 
the same way as before. 


dcr 

dJ 


Cf . — 


— Qfo 


dz 

yl-2e 


ii 1 2 ,coll J 2 ,finite 

cl/ dcr A dcr A dcr A 


(80) 


where I is given by Eq. m- One sees immediately that the soft contribution is exactly the 
same as for the (/) term: 


_ Uf _ _U'F 


dJ 


dJ 


(81) 


As to the collinear piece, we note that in the collinear limit Z = 0 the jet distribution reduces 
(by applying Eq. (|49d[) ) to Sf\ and one getsQ 
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zVgq{z, e)S} (fc; x{l - z)) - 2S) {k; x) 


P) I dx fq^\x) X 
0 

(2)7 


^The UV cutoff A for the initial state collinear singularity is in principle independent of the one adopted in 
the final state collinear term of Eq. ([t^. We use the same cutoff for all collinear subtractions in view of its 
identihcation with the factorization scale A = fj,F- 

























With the change of variable z s- 1 — z, and performing the I integration, we obtain 


J 2,coll 

dJ 


= I dk hf\k)hf\k) dx /f (x) X 

I [1 + 2 :' + £(1 - zf] (fc; xz) , 


(83) 


where the standard ()+ regnlarization has been nsed. To separate singnlar from hnite pieces one 
has to perform an ^-expansion inside the ^-integral. In order to cast the collinear singnlarity 
into the standard form, we introdnce the fnll 4-dimensional Altarelli-Parisi q ^ q splitting 
function^] 
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We then can write 
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In the last expression the hrst term is again the LO jet cross section, mnltiplied by a singnlar 
factor; the second term contains the proper qnark corrections to the qnark distribntion fnnction, 
while the remaining pieces are hnite in the e —0 limit. 

The last contribntion in the Cf part is regnlar in 4 dimensions and dehnes another hnite 
term 

1 i,finite A o- P AI 

jf- -y ^[/(.-,!)-/(0,()-(/(.-,0)-/(0,0))e(A'^-C)] (86) 

and will be nsed in onr hnal resnlt. 


4.4.2 CA term 

The term proportional to Ca in Eq- O reads 
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It shows a g II 1 collinear singnlarity corresponding to the k' = d pole. Becanse of the nnmerator, 
the fchintegration is not really singular at q = 0 or at qr — zfc = 0, except for z — > 0. In fact, in 

®The appearance of a g —> q splitting function is simply due to the relation ^^^(l — z) = Pqq{z) valid for 
z < 1. 
















the high energy limit (j^), q hxed and z -C 1 correspond to gluon 3 being in the central region, 
where we have 


g- i.Q-zk) 
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This is exactly the expression entering the differential partonic cross section (^2]) in the central 
region that provides the LL contribution. In other words: Eq. (^7|) contains the O' || 1 collinear 
singularity in the whole z-range; for finite values of z neither the (? || 3 nor the 1 || 3 collinear 
singularities are really present. However, in the (gluon) central region z -C 1 their “collinear 
denominators” degenerate, providing the soft singular real part of the leading kernel l/q\ 
The jet distribution functions will become essential in disentangling the collinear singulari¬ 
ties, the soft singularities, and the leading log s pieces. The basic mechanism can be understood 
as follows: 


• When the outgoing quark 1 is in the collinear region of the incoming q, i.e., yi —> oo, 
quark 1 cannot enter the jet; only gluon 3 can thus be the jet, 1/3 is fixed and no logarithm 
of the energy can arise due to the lack of evolution in the gluon rapidity. No other singular 
configuration is found for the quark when J = {3}. 

• In the composite jet configuration, i.e., J = {1, 3}, the gluon rapidity is bounded within a 
small range of values, and also in this case no logs can arise. There could be a singularity 
for vanishing gluon 3 momentum: even if the 1 || 3 collinear singularity is absent, we 
have seen that, at very low z, a soft singular integrand arises. However, the divergence is 
prevented by the jet cone boundary, which causes a shrinkage of the domain of integration 
~ for z ^ 0 and thus compensates the growth of the integrand. 


• The one-quark jet configuration J = {1} allows the gluon to span the whole phase space, 
apart, of course, from the jet region itself. The LL term arises from gluon configurations 
in the central region. But also here, like in the negative rapidity region discussed in 
Sec. it is crucial to understand to what extent the differential cross section provides 
a leading contribution. It turns out that the coherence of QCD radiation suppresses 
the emission probability for gluon 3 rapidity j /3 being larger than the rapidity yi of the 
outgoing quark 1, and an angular ordering prescription similar to that of Eq. (^) holds. 
This will provide the final form of the leading term, i.e., the appropriate scale of the 
energy and, as a consequence, a finite and definite expression for the one-loop jet vertex 
correction. 


As a first step we isolate, in Eq. (^), the initial state O' || 1 collinear singular contribution 
and define, like in the C^^-term analysis, the collinear term by setting A;' = 0 (except in the 
pole), and by introducing an UV cutoff: 
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In this expression the jet distribution, because of Eq. (|49c|) , reduces to {k; xz). By including 
the Me A constant factors and the l/k‘^ pole, we reconstruct the gluon-initiated LO jet vertex 
(see Eqs. (^) and (|I9[)). Note also that 
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where Pgq{z) is the 4-dimensional q ^ g splitting function. It will be used to dehne the quark 
correction to the gluon distribution function, and the lower bound Zcut can be set equal to zero 
(up to an error of the order t/s). The transverse integral is easily performed, and we obtain 

'A 2 ' 
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TT 
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By subtracting the collinear term (^) from the total Ca term (^), we obtain an expression, 
that is hnite in the e —> 0 limit, and whose value increases logarithmically with s because of 
the dz/z integration with lower bound Zcut ~ s~^P. As already mentioned, it is extremely 
important to understand the extension of the phase space contributing to the cross section at 
NLL accuracy, i.e., by including also the constant terms. That is actually crucial in the small-z 
region where the integrand develops a soft singularity which, if not properly absorbed in the 
real part of the kernel, could give rise to spurious divergencies spoiling the whole procedure. 

It is instructive to review how the mechanism of suppression of the differential cross section 
sets in when the gluon is emitted at an angle smaller than that of the quark. At hxed quark 
momentum pi, i.e., at hxed transverse momentum k' ^ 0 and longitudinal momentum pf ~ 
(1 — z)xy/s/2 — one can imagine J = { 1 } hxed by the jet condition —, and at hxed gluon 
transverse energy = \q\ and longitudinal momentum ~ xz-\/s/2, we perform an azimuthal 
average of the subtracted Ca diherential cross section with respect to the angle 03 of the gluon. 
The collinear subtraction actually does not contribute, because it refers to conhguration with 
k' = 0. If in this averaging procedure we neglect the variation of hf\k), it is sufficient to 
calculate the azimuthal average of the factor 


{l-z)q- {q- zk) 


= - 0(^3 - Z{E^ + E^)) 


q^{q-zkf r 

This relation is exact and clearly shows that, outside the angular ordered region 
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(93) 


2 ; 1 — z 

there is no contribution to the cross section. In practice, by taking into account the variation of 
h® (fc) during the averaging procedure, we have to replace “no contribution” by “suppression”. 
Note however that, in the limit q ^ 0 (which includes the soft region), the variation of k' goes 
to zero as well, so that Eq. (|^) is really an accurate statement in the “dangerous” part of the 
phase space. Moreover, Eq. (^) shows that the 1/q^ kinematic dependence of the LL kernel 
governs the differential cross section up to the very end of the angular boundary. Therefore, 
we dehne the LL contribution in the “lower half region” i/g > 0 by 
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where we have imposed the J = {1} jet condition. 

The remaining part of the Ca term is hnite in 4 dimensions and constant in energy, so that 
we can set e = 0 and Zcnt = 0 to dehne the constant part 

dj ■ [dj dj dj J 


(95) 












5 The NLO jet vertex: sum of real and virtual correc¬ 
tions 


Having completed the calculation of both the virtual and real corrections in the whole phase 
space, we are going to collect all partial results and to show that the complete one-loop jet cross 
section can naturally be htted to the form of Eq. (^). Table ^ summarizes the decomposition 
of the one-loop jet cross section and gives the references of the various contributions. 

In Sec. 1^ we have presented the virtual contributions to the jet cross section which, after 
renormalization of the coupling, assume the form of Eq. (|4^) . We have already taken into 
account the contribution coming from the H;, impact factor correction in Sec. |4.3| by combining 
it with the “upper half region” real contribution in Eq. (|68|) . The remaining virtual terms can 
be conveniently rewritten in the form 
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(96a) 
X (96b) 


where the virtual kernel has been dehned in Eq. (p9| ) and the coefficient Hg is given in 

Eq. (^51) . The energy scale sq in Eq. ( |96a|) is constrained to satisfy so(fc, k) = k^ = Ej, because 
the delta function inside the virtual kernel sets k = k', and the delta function in the LO jet 
vertex sets k'"^ = Ej. The general form of the energy scale for k' ^ k will be hxed in a moment 
by the real LL contribution. 

they show the same structure and 


We hrst join the LL real contributions (|68a|) and 
differ only in the logarithmic term. The sum of the two logarithms yields 
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The denominator in the argument of the logarithm dehnes the energy scale sq- However, 
as we have already pointed out, there is some freedom in choosing the LL subtraction and, 
correspondingly, the denominators in the logs. We can obtain a more symmetric expression 
by dehning the LL term in the “upper half region” as suggested at the end of Sec. by 


Table 1: Schematics of the decomposition of real and virtual one-loop corrections to qb scat¬ 
tering and references of the corresponding equations. 
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replacing the scale max(£' 2 , E^) by E 2 + -E 3 . This amounts to using the same prescription for 
defining the LL contribution in the “upper half region” 1/3 < 0 and in the “lower half region” 
1/3 < 0. With this choice, the full LL contribution to the jet cross section, including the virtual 
correction 
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( 100 c) 

We identify Eq. 


These constraints are consequences of the QCD coherence effects 
with the first term of Eq. (^^. 

Let us stress that the energy scale in Eq. (p9| ) — and in general all the scales satisfying 
Eqs. (|100|) — arises naturally when one requires impact factors and PDFs to have standard 
collinear properties and the remaining non-leading-log term (p5|) to be hnite in both the physical 
£ —>■ 0 and high-energy s ^ 00 limits. Choosing a scale of the energy outside the class dehned by 
Eqs. (p.00|) while preserving the above properties, requires the introduction of additional NLL 


operators (see Eq. (|71|)) which has to be added as multiplicative corrections to the Green’s 
function. If, for instance, we adopt Sq = |fc||fc^|, then the Green’s function in Eq. ( p3b|) has to 
be replaced by 
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The second term of Eq. (P^f ) has already been obtained in Eq. ( |t) 8 b| ). We remark again that, 
with the choice of the energy scale (PUj), the one-loop impact factor correction is no longer 
given by Eq. (|69|) but differs by a finite part. However, its actual expression is irrelevant for 
the jet vertex. 

We now consider the sum of all e-divergent contributions that have not been included in 
the impact factor term ( | 68 b|) . They can be found in Eqs. (|96b| j77|j7^JM|j85| J^) and add up 
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All double poles have cancelled out. Single poles only appear in connection with the splitting 
functions: they are shown in the hrst two lines of Eq. ( |103|) , and they contribute to the third 
term of Eq. (^) which contains the PDF one-loop corrections. Note however that, since the 
present analysis has been restricted to the case of incoming quarks, we have obtained only the 
quark-initiated corrections to the quark and gluon PDF, i.e., only the term c = g in the sum of 
Eq. (pbl). The gluon-initiated corrections require an incoming gluon out of hadron H and will 
be presented in a forthcoming paper [^. It is also clear that the cutoff A can be identihed 
with the factorization scale /j^f- 

The last two terms of Eq. (|103|) are regular at £ = 0 and can be combined with the hnite 
parts (17^), (jS^) and (pSl). The resulting expression can be cast into the form 
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which dehnes the NLO correction to the quark-initiated jet vertex 
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It clearly depends on the jet dehnition and on three scales: the energy scale Sq (via the sub¬ 
traction of the LL term cx 1/^), the factorization scale A = /ij? and the renormalization scale 

/i. 

Eqs. (|104|) and (|105|) provide the fourth and last term of Eq. (|35D and represent the main 
result of our study. 


6 Conclusive remarks 

In this paper we have investigated a particular class of jet hnal states in the high energy region. 
Both the theoretical and phenomenological motivation comes from the interest in the Regge 
limit of QCD: the jet production processes - forward jets in ep collisions and Mueller-Navelet 
jets in hadron-hadron collisions - have been designed to investigate BFKL dynamics at present 
and future colliders. Whereas the BFKL Pomeron is known at LG and NLO, a complete NLO 
analysis and comparison with data has not been possible, since neither the photon impact factor 
nor the jet vertex have been calculated in NLO. On the other hand, there is no doubt that at 
moderately large energies, NLO corrections to the asymptotic LO behaviour are important for 





























a reliable description of the processes under consideration. It is the purpose of this paper (and 
a forthcoming one) to calculate the jet vertex in NLO: in this first part we have studied the 
quark-initiated jet vertex, whereas the gluonic counterpart will be presented in a subsequent 
paper. 

In order to extract the NLO jet vertex we have calculated the cross section of the process 
quark b+quark a —>• quark b+X+jet at order apart from the NLO corrections to the impact 
factor of quark b and the contribution to gluon production in the central region, this process 
provides the NLO corrections to the quark-initiated vertex quark a jet. As an important 
theoretical result we have verified that the factorization form (p3a|) holds: this proof follows 


from the fact that we have been able to separate, in the sum of virtual and real corrections, the 
collinear singularities which go into the renormalization of the parton density of the incoming 
quark a, and the gluon emission in the central region which is part of the LO calculation. 

Another theoretical issue of interest is the dependence upon the scales. In leading order, 
results of the BFKL calculations and the jet vertex are insensitive to both the energy scale Sq 
and to the renormalization scale /r. It is only at the NLO level that the dependence on these 
scales is being determined. 

The central result of our calculation is the expression ( p.05| ) for the NLO jet vertex. Using 
the factorization property ( P3a|) , we can use our result also for the upper incoming quark b, i.e., 
for the ‘symmetric’ Mueller-Navelet jet production process q + q ^ jet + X + jet. As the final 
step, we will have to allow for the production of an arbitrary number of gluons between the jets 
which is described by the NLO BFKL Pomeron. This step will be presented in the companion 
paper containing the NLO gluon-initiated jet vertex. 

With our expression for the jet vertex we have provided a finite integral that can be com¬ 
puted numerically, e.g. via Monte-Carlo integration. Such a numerical study requires to specify 
the jet algorithm described by the function 5®: so far we been rather general, but clearly the 
numerical results will depend on the choice of the jet algorithm. 
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